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Abstract. We consider inverse problems for wave, heat and Schrodinger- 
type operators and corresponding spectral problems on domains of R n and 
compact manifolds. Also, we study inverse problems where coefficients of 
partial differential operator have to be found when one knows how much 
energy it is required to force the solution to have given boundary values, 
i.e., one knows how much energy is needed to make given measurements. 
The main result of the paper is to show that all these problems are shown 
to be equivalent. 



1 Introduction 

The goal of this paper is to explain the equivalence of various types of data 
used in inverse boundary value problems. Our main interest lies in the in- 
verse boundary spectral problems and the inverse boundary value problems 
in the time-domain for equations with time-independent coefficients. There 
are numerous examples when a solution of a particular inverse problem is 
used to solve an inverse problem of another type. For example, A. Nachman, 



J. Sylvester, and G. Uhlmann [NSU] solved the inverse boundary spectral 



problem for a Schrodinger operator by reducing it to the inverse boundary 
value problem in fixed frequency and then using the method of complex geo- 



metric optics [3U|. Similarly, inverse problems in the time-domain are often 
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reduced to problems in the frequency domain (see e.g. Q). In this paper 
we are interested in the equivalence of different types of boundary data for 
inverse problems and describe some procedures for transforming between 
these problems. We believe that the equivalence of various inverse problems 
is useful from the theoretical point of view since it enables us to translate the 
results obtained for one inverse problem to other problems. The equivalence 
is also useful for applications as it makes possible to use reconstruction al- 
gorithms developed for some particular type of inverse problems to different 
inverse problems. Moreover, in various applications some measurements are 
considerably more difficult to make when compared to other measurements 
which, however, can give an equivalent information. For instance, it is of- 
ten difficult to measure the phase of a physical field but the energy which 
is required to force the boundary value of the field to a given one is often 
known very precisely. In theoretical inverse problems this idea goes back to 
A. Calderon who in his seminal paper of 1980 considered an inverse problem 
for the conductivity equation from the point view of energy measurements. 
Similarly, the energy or interference based measurements have been used in 



many applications, e.g. in impedance tomography (see e.g. [GIN]) and near 
field optical tomography (see e.g. [ |SM| |). 

We start with an elliptic problem in the frequency domain. Let f2 C R' m 
be a smooth bounded domain and a(x, D) be an elliptic 2nd-order partial 
differential operator, 

a(x, D)v{x) = -g-^dj (g l/2 g lk (x)d k v(x)) + q(x)v(x), (1) 

where [g^ k (x)] is a positive definite smooth real matrix in Q, [gij(x)] is 
the inverse matrix of [g^ k {x)\, g{x) = det(gij), and q(x) is a smooth real 
function. As usual we use Einstein's summation over repeated upper and 
lower indices. The operator a(x, D) is a Schrodinger operator in f2 which 
corresponds to the metric tensor gij. 

Two objects related to the Dirichlet problem are the boundary spectral 
data and the Dirichlet-to-Neumann map A z . Namely, if A is the operator 

Av = a(x, D)v 

with 

V(A) = H 2 (n) PiH^n), (2) 
then its boundary spectral data is the collection 

{A,, Btpilea : 1 = 1,2,...}, (3) 
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where A; and <pi are the eigenvalues and the normalized eigenfunctions of 
the operator A, 

a(x, D)tpi = Xiipi, w\ 9n = 0. 

The boundary operator B is given by 

B<pi = gWcpVjdwu (4) 

where v = (ui, . . . , u m ) is the interior unit normal vector to dO, in the Eu- 
clidean metric. To clarify the meaning of the boundary operator B, we note 
that if n is the interior unit normal vector to d£l with respect to the metric 
g, then Bu = pd n u, where the weight p is given by 

p(x) = det^W^V^)^) 172 - (5) 

Here ggn is the restriction of the metric tensor gjk(x) to the tangent plane 
H = T x (dQ) of dQ at x. Thus, integration by parts gives rise to the formula 



/ 



a(x,D)vudV g = / (g lj djV d k u + qvu) dV g + / Bv\gn u\gndS. 
Jn Jan 



Here dS is the Euclidean area element on d^l while dV g = g 1 l 2 dx is the 
Riemannian volume element on in the metric g. 

Then, for 2^ A;, the (fixed-frequency) Dirichlet-to-Neumann map A z is 
defined by 

A z h = Bv h z \ dn , h£C°°(dn), (6) 

where v\ solves the boundary value problem 

a(x,D)v h z =zv h z , v% n = h. (7) 

We are interested in the case when A z is given at all possible values of the 
spectral parameter z £ C. This data is called the Gel'fand spectral data. 

Operator a(x, D) defines in a natural way initial boundary value prob- 
lems for the wave equation, 

(<9 t 2 + a(x,D))u f = in !JxR + , (8) 
n / |anxR + = / € C^(dQ x R+), u f \ t=0 = d t u f \ t=0 = 0; 

for the heat equation, 

(d t + a(x,D))w f = in n x R+, (9) 
w f \dnxn+ = f, w f \ t =o = 0; 
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and for the non-stationary Schrodinger equation, 

(id t + a(x,D))ip f = in O x R + , (10) 
ip f \dnxR+ = f, ip f \t=o = 0. 

For the initial boundary value problems (p|)-(|i"o|) we define the non-stationary 
Dirichlet-to-Neumann maps (response operators) R w , R , and R s which 
correspond to the wave, heat, and non-stationary Schrodinger equations, 

R w f = Buf\ dQxR+ , (11) 
R h f = Bwf\ dnxR+ , (12) 
R s f = B?P f \ dnxR+ . (13) 

We note that R w and R h are real-valued operators, i.e. R w /, R h f are real 
as soon as / is real. Later, this will allow us to restrict ourselves to the 
consideration of only real f's when dealing with problems (||) and (||) . 
There is a natural concept of energy for the wave equation given by 

E w {u,t) = ^ f (^ fc a^(x,t)a fc ^M) + '7(^|n(x,t)| 2 + |a t n(x,t)| 2 ) dV g (x). 

Similarly, for the non-stationary Schrodinger equation we define the energy 

E s &,t) = ^J^ k djij(x,t)d k W^)+Q^)\^,t)\ 2 ) dV g (x). 

The notion of energy (or, better say, heat content) for the heat equation has 
sense in the case q = and is given by 

E h (w,t)= [ w(x,t)dV g (x), (14) 

with the term qw physically corresponding to a heat source. For the wave 
and non-stationary Schrodinger equations the energy is preserved in the 
absence of the boundary and interior sources. In the case of the initial 
boundary value problems with / 6 C§°(d£l x [0, T]) this means that the 
energies E w (u f ,t) and E s (ipf , t) are constant for t > T. Therefore, the 
energies E w (u-f ,t) and E s (ipf,t) are brought into Q through the boundary 
d£l x R + . We define the total energy fluxes through the boundary, U w (f) 
and IT S (/) for the wave and non-stationary Schrodinger equations as 

U w (f) = lim £<V,i), IF(/) = ton E s (^,t), (15) 
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where / G Cfi°(dfl x R+). Actually, IL w (f) and II s (/) are given by some 



quadratic forms of / (see formulae (28) and (p5|)). The total energy fluxes 



are important concepts, as they are often easier to measure in practice than 
the corresponding non-stationary Dirichlet-to-Neumann maps. 

Now we are in the position to formulate various inverse problems in the 
domain 0, which are related to the above concepts. 

Elliptic problems: 

i. Given boundary spectral data, {A;, Bifi\g^ : I = 1, 2, . . . }, determine 
gi k and q. 

ii. Given Gel'fand spectral data, {A 2 , z S C}, determine gi k and q. 

Hyperbolic problems: 

hi. Given a hyperbolic Dirichlet-to-Neumann map R w , determine gi k and 
Q- 

iv. Given a hyperbolic energy flux TP", determine g^ k and q. 
Parabolic problems: 

v. Given a parabolic Dirichlet-to-Neumann map R h , determine gi and 

q- 

Non-stationary Schrodinger problems: 

vi. Given a Dirichlet-to-Neumann map for a non-stationary Schrodinger 
equation, R s determine gi k and q. 

vii. Given an energy flux II s for a Schrodinger equation, determine gi k and 
Our main result is: 

Theorem 1 Inverse problems i.-vi. are equivalent, i.e. any of the data i.- 
vi. determine all other data. In the case of problem vii. we should require, 
in addition, that A = is not an eigenvalue of A. Then vii. is equivalent to 
i. — vi. 

We note that if S cr(A), where a (A) = {X\ < A2 < • • • } is the spectrum 
of A, then II s determines all {A;, Bipi\gn} for A/ 7^ 0. Similarly, II s in this 
case determines all Dirichlet-to-Neumann maps (||), (JTT|) -(13) upto a finite- 
dimensional operator. 
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The authors want to emphasize that some of the above equivalences are 
known at least on the formal level. Rigorous proofs of some of them are given 
in [KKL|, [Q. In these cases, we will describe the corresponding results 



very briefly and refer to KKL |. However, in particular for the Dirichlet 
problem considered in this paper, care should be taken to make these formal 
constructions rigorous. 

Although solving inverse problems i.-vii. is not the main goal of this 
paper, we would like to mention the following application of Theorem [l|: 

Theorem 2 Assume that and any of the boundary data i.-vii. are given. 
Then it is possible to construct gi k and q upto a diffeomorphism, that is, we 
can find the equivalence class 

[a(x, D)] = {&*(a(x,D)) : <I> : f2 — > is a diffeomorphism, &\dn = id\oti}, 
where <$* is the pull-back of the diffeomorphism <3?. 

For problem i. the proof of Theorem can be found in [Kl], [K2|, and, 



in more detail, in [KKL|, Section 4.5. Moreover, it is shown in [KK] that 
Theorem ^ is valid for the case when we know all but a finite number of 
boundary spectral data. Other cases follow then from Theorem |] and the 
remark after it. 

In Sections 2-4 we will describe some constructive methods to transform 
data i.-vii. into each other. These methods require neither solving the cor- 
responding inverse problem nor using an analytic continuation. We will not 



provide all technical details by either referring to [KKL] or Appendix 1. In 
Section 5 we will discuss some generalizations of problems i.-vii. In partic- 
ular, we consider the case when the data is given on a finite time interval 
or a(x, D) is different from ([!]) with some technical details concerning the 
energy flux for the wave equation considered in Appendix 2. 

2 Equivalence of the boundary spectral and Gel'fand 
data. 

i — > ii . Given {Aj, B(pi\gQ, I = 1,2,...} one can formally construct the 
Dirichlet-to-Neumann map A z , 

(A^)(x) = f^^M f h {y)B^(y)dS y , (16) 

formally ^ Z — A\ Jg n 



1=1 
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where dS y is the area element of the surface d£l in R m and, without loss of 
generality, we take real <pi. However, the right hand side of (16) does not 
converge [NL|. Therefore, the sum in © has to be regularized. This was 
done by A. Nachman, J. Sylvester and G. Uhlmann in NSU | for the case 
gi k = S^ k , i.e. a(x, D) = —A + q. Their construction was based on the fact 
that — A + q is a relatively compact perturbation of the Laplace operator 
—A and the boundary spectral data and the Dirichlet-to-Neumann maps for 
the Laplace operator in a given domain Q can be always found. However, 
—g 1 ^ 2 dj(g 1 ^ 2 g^ k dk) +q with unknown gi k is not a relatively compact pertur- 
bation of some known operator and thus one needs a different regularization. 
To this end, we start with the eigenfunction expansion 

oo ^ 

^ = J2^^ ^ = — -T-(h,B<pi) L 2 {m) , (17) 
1=1 1 

which converges in L 2 (Q). Differentiating equation ( |l7| ) with respect to z 
we obtain a formula for the derivative d z v z , 

OO - 

d z v z=-y2l rT2(KBtpi) L 2 {m) tpi. (18) 

1=1 ^ Z ~ XI > 



Now the sum in the right hand side of (|i~8| ) converges in H 2 (Q). Thus, 
the boundary spectral data determines the derivative of the Dirichlet-to- 
Neumann map, 

oo ^ 

L 2 (dn)B(pi\dn, (19) 

where the right hand side converges in H x / 2 {d9). However, A z h can be 
represented as 

A z h= lim (A_ r2 /i+ / d z >(A z ,h) dz' ) . (20) 
r^oo y J_ T 2 J 



Thus, (|19|), ( po| ) determine the map A^ as soon as we can find the asymp- 
totics of A_ r 2 h, t — » oo. It is shown in Appendix 1 that 

(A_ T 2 h){x) = -T P {x)h(x) - -p(x)H(x)h(x) + 0(r -1 ), r -> oo, (21) 

where p is given by formula (|5|) and H (x) is the mean curvature of d£l with 
respect to the metric g^. In turn, the mean curvature can be found from 
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the symbol of d z A z which is determined from the boundary spectral data by 



( |i9|) (for a detailed construction see Appendix 1). Thus formulae (|19|)-(21) 
determine A z from the boundary spectral data. 

ii. — > L Formulae (|l8|)-(|l9|) show that A z is a meromorphic operator- 
valued function which has simple poles at the eigenvalues z = A/. This 
implies that A z determines the eigenvalues A/. Moreover, one can show that 
the residue of A z at z = A/ is a finite dimensional integral operator with the 
kernel 

Kfa, y) = J2 B^ k {x)B^ k {y), x,ye cKl (22) 
keLi 

Here Li is the set of integers k such that A& = A;. Clearly, when A; has mul- 
tiplicity one, the kernel Ki(x,y) determines Bifi{x) up to a multiplication 
by ±1. In general, the linear independency of Bip^x) makes it possible to 
find B(fk(x), k € L[ up to an orthogonal transformation. More precisely, we 
can find 

tj( x ) = ^2 a jkBip k (x), j e L h (23) 
keLi 

where [ctjk] is an orthogonal matrix (for details see KKL|| ). 

Clearly, £j are themselves the boundary values of some eigenfunctions 
tpj which form another orthonormal basis of the eigenspace corresponding 
to the eigenvalue z = Xi, i.e. £j = Biftj . Thus, the Gel'fand data determines 
the boundary spectral data. 

3 Equivalence of the hyperbolic and elliptic bound- 
ary data 

ii — > Hi. To construct R w from A z , z £ C, let 

f(x,u)= / e-^f(x,t)dt (24) 
J o 

be the Laplace transform of / G Cq°(3J1 x R+). It is well defined for to G C 
and satisfies 

\\J{u)\\ H z/2 {dn) <C N {l + \u\)- N , Re^>0, (25) 

for any N > 0. Let v(uj) be the solution of the Dirichlet problem (^) with 
z = —u! 2 and h = f(ui). Since 

\\v h z\\HHn) < C(l + \z\) (1 + dist- 1 (z,a(A)))\\h\\ (26) 



S 



estimate (25) implies that v{uj) decreases rapidly when \uj\ — > oo in the right 
half-plane, Reuj > \x > 0. Moreover, Bv(u)\qq = A_ w 2/(o;). 
Using the inverse Laplace transform, we see that 



1 



oo 



u (x, t) = ^- I e^+^vCx, \i + is)ds 



DC 



is well defined in C°°(R + , H 2 (Q)) for fi > y/JIo, where 

Ho = max{0, -Ai}, (27) 
and solves initial boundary value problem (jg). Therefore, 

1 1 00 ~ 
#7 = Su| 9 nxR + = ^ y A_ (/i+is)2 /(x, /x + is)ds. 

Thus A z determines the map R w . 

Hi. — > iv. It follows from the definition of U w that 

/» poo 

U w f=Re / (R w f)(x,t)d t f(x,t)dtdS x , (28) 

which implies that the hyperbolic Dirichlet-to-Neumann map determines 
IT". 

iv. — > i. The quadratic form U w (f) determines the bilinear form 

IT [f,h] = ^Re y y°° ((R w f)(x, t)d t h(x, t) + d t f(x, t)(R™h)(x, t)) dtdS^. 

Because i?™ maps the real valued functions to the real valued functions, we 
can separate the real and imaginary parts of H w and construct the complex 
bilinear form 

[/,/»] = J Re / [°° ((R w f)(x,t)d t h(x,t) + d t f(x,t)(R w h)(x,t)) dtdS x . (29) 
1 JdQ JO 

Using the Plancherel formula we can represent the integral in the right-hand 
side of ( |29| ) in terms of the Fourier transforms /, h of /, h, where e.g. 

f(x,k) = [ e ikt f(x,t)dt. 

Obviously, R w f(k) = A k 2f(k). As A k 2 has poles at points k = ±a/A7, the 
contour of integration in the Plancherel formula is shifted in the complex 
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Figure 1: Contour L 
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domain to the boundary of the strip {k £ C : \Jm k\ < fj,} (see Fig. |l|) 
and we obtain 

mf,h] = ^- [ [ A k2 f(-k)h(k)dS x kdk (30) 
47 ™ Jr M Jan 

(for details see [ KKL |). 

Taking into account the analyticity of f(k) and h(k) and the repre- 
sentation ( p2| ) of the residues of A 2 at z = A;, estimate (p6[ ) and Weyl's 
asymptotics for eigenvalues, the residue theorem yields that 

+ (/(-A)»^i) (fcCAO.B^-)) . (31) 



This formula provides a possibility to find the kernels Ki(x,y) given by 
formula (|22|). Indeed, let / = f T , h be of the form 

f T (x, t) = F(x) X (t - r), h(x, t) = H(x) X (t), (32) 



where F,H& C°°(dn) and \ G C^°(R + ) is real-valued. Then, formula @ 
yields that 

1 °° 

ng[/r s /i] = -^cos(yA7r)|x(yA7)| 2 (^F,^) L2( ^), (33) 

Z=l 
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where K\ are the integral operators with kernels Ki(x,y). Analyzing the 
asymptotic behavior of the sum (|33| ) when r — ► oo, we can find A/ and 
(KiF,H). As F and H are arbitrary, we find the kernels Ki(x,y) for any I 
(for details see [|KKL| ) . Due to the implication %%. — ► i., this also determines 
the boundary spectral data. 

4 Equivalence of the parabolic and non- stationary 
Schrodinger boundary data and the boundary 
spectral data 

ii — > vi. Let / G C^°(9f2 x R+). Its Fourier transform with respect to time 
variable t, f(x, k) is a rapidly decreasing function in the upper half plane 
C+ = {keC: Imk > 0}. Let 

rj(x,t) = — [ e~ lkt rj(x,k) dk, fi>0, 

where 

a(x,D)i](x,k) = -krj(x,k) in fi, r)(x, k)\gu = f(x, k). 

Using estimate ( pq) , we see that rj(x,t) is the solution of the initial boundary 
value problem (|1C|), i.e., r](x,t) = ipf(x,t). As Brj{k) = K^f(k), we obtain 

R S f = ^l e- lkt A^f(k)dk. (34) 



vi — >■ vii. By direct calculations we obtain 

p roc 

n s (/)=Re/ / R s f(x,t)d t f(x,t)dtdS x . 
Jan Jo 

Thus, R s determines II s . 

For the future we note that R s determines also 

||^(oo)|U 2 = lim U f (t)\\ L2 . 

t— >oo 



Indeed 

II,/,/ 



(35) 



(oo)\\ L 2=i[ [ (R s f(x,t)f(x,t)-f(x,t)R s f(x,t))dtdS 
Jan Jo J 
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P f'QO 

2 1m / / f(x,t)R s f(x,t)dtdS. (36) 
Jan Jo 

This implies that R s determines also the energy flux n|; o (/) which corre- 
sponds to the potential q + Eq, Eq £ H. 

vii — > i. Quadratic form II s (/) determines the form Uf,[f,h], 

U s r [f, h) = ~ Re f [ (R s f{x, t)dth(x, t) + d t f(x, t) R s h(x, tj) dtdS x = 
1 Jan Jo v ' 

^Re I f (R?f(x, t)dth(x, t) + d t f(x, t) R s h(x, t)) dtdS x . 
Taking / = if in this formula we obtain also that 

KilM =~lm f (R s f(x, t)dth(x, t) + d t f(x, t) R s h(x, t)) dtdS x 
Thus, the quadratic form II s (/) determines the bilinear form II s [/, h], 
n s [/,/i] = \f (R s f(x, t)dth(x, t) + d t f(x, t) R°h(x, i)) dtdS x (37) 

= \ I [ (e-^ R s f(x,t) e^d t h{x,t) + e^d t f{x,t) e-^R s h(x,t)) dtdS x , 
2 Jan Jo ^ ' 

where \i > is arbitrary. Then, by the Parseval identity, 

n s [M] = \! I" ((e-^R s f) A (x, k)(el*dth)*(x, k) + 
1 Jan J-od v 



+ (e'" t d t f) A (x, k)(e-^R s h) A (x, k)\ dkdS x . 
Using formula (p3|), we see that 

n s [/, h] = \r \-i(k - in){f(k - ifj), {k^fi){k + in)) 

^ J -00 

+i(k + iy) {(^-k~i^7){k + ifj,),h{k - i/i))| dk = 

= ^-[ [ zA_ z f(x,z)h(x,z)dzdS x , (38) 
4? ™ Jan Jr„ 
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where is given in section ||. In deriving the last equation in (38) we use 
the identity (A_^)* = A_^. 

Since A_ 2 has poles at points z = —A;, we can apply the residue theorem 
as in section [3| and obtain 



1 oo 

n s [f,h] = - Y,MKif(-Xi)M-^i))- 

1=1 

Taking / and h of form (||) with F,G e C°°(dn), we get 

^ oo 

U s [f,h] = -Y,^ieM-^ir)(KiF,H)\x(Xi)\ 2 . (39) 
1=1 

The sum ( |39"| ) considered as a function of t £ R determines all the eigen- 
values A/, A; ^ and the corresponding bilinear forms K\ (see section |for 
further details). Thus II s determines the boundary spectral data except for 
Bifi with A/ = 0. 

We note that the terms Bip\ for which A/ = do not appear in series 
( |39| ) and thus the boundary spectral data corresponding to the eigenvalue 
can not be determined from the asymptotics. 



vi. — » i. As noted after formula (|3q), having R s it is possible to construct 
the form 11^ . Thus, by choosing sufficiently large Eq and using formulae 
(|35|), (|36| ) we can determine all the eigenvalues A/ and corresponding func- 
tions Bipi. Similarly if, in addition to the energy flux II s , we know the 
quadratic form / i— ► \\i/)f (oo) || 2 , we can also determine all A;, Bipi. 

ii. — > v. Using the Laplace transform instead of the Fourier transform 
and following the same steps as for the non-stationary Schrodinger operator, 
we obtain 

R h f( x ,t) = ^- [ e^J(u)dLJ. (40) 

(compare with formula (j34|)). 

v. — > i. Similar considerations as for the non-stationary Schrodinger op- 
erator show that R h f(x,uj) = A_ a ,/(a;, uS) for Rew > fio where fiQ is of form 
(p7|). Using, as earlier, the boundary sources of form ( |32|) and applying the 
Parseval formula and the residue theorem, we obtain 



I dn Jo 



poo 00 

/ f(x,t)R h h(x,t)dS x dt = y2e- x < T \x(\ l )\ 2 (K l F,H) (41) 
Jo i=i 
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(for details, see [ |KKL| ] ) . As above, the asymptotics of the sum in the right 
hand side of (f4l|) for r — » oo determines the boundary spectral data. This 
proves Theorem [j]. 



5 Generalizations 
5.1 Continuation of data. 

From the point of view of uniqueness, we are often interested in the case 
when the data is given only on a finite time or spectral intervals. In the 
other words, we are given only the restrictions R w f\dQx[o,T], R s /|dfix[o,T] 
or R h f\dnx[o,T] on the time interval [0, T] for the time domain problems 
or A z for z £ I where I <Z C for the spectral problem. In particular, as 
inverse problems are usually ill-posed, it is important to find out whether 
these restricted data determine the complete data, i.e. R w , R s , R h or A z 
for z G C. 

In the spectral problem we see that A z is a meromorphic operator- valued 
function of z. Thus, by analytic continuation A z , z £ I determines A z for 
all z G C as soon as / has an accumulation point. For the time-domain 
problems, we can use a complexification of time for the similar purposes. 
Indeed, when / G C£j°(dtt x [0,To]), we can use the Paley- Wiener theorem 
and estimate ( |26| ) to deform the countours of integration in integrals ( |34| ) 
and ( [40D with t > Tq to the contour T^, \i > hq. This contour consists of 
a straight ray (— oo — ijx, —in), semicircle of the radius \x and straight ray 
(ifjL, — oo + ifx). Then for any function / G C^°(dQ x [0, To]), the functions 

t^R h f(x,t), t£{Ret>T }, 
t^R s f(x,t), t G {Ret > T , Imt > 0} 

are holomorphic functions. Moreover, R s f(x,t) is continuous with respect 
to t G C up to the real semiaxis t > T. Thus, for any / G C^°(dn x [0,T ]) 
R h f\dQx[o,T] an d R s f\dQ X [o,T]i T > T determine, correspondingly, R h f and 
R s f on dU x R + . However, R s and R h are linear operator which commute 
with the time delay operator Y T , Y T f(t) = f(t — r). Thus, the restrictions 
R h f\dnx[o,T] and R s f\ d nx[o,T] with / G Cg>(dn X [0, T ]) determine R h f and 
R s f for any / G C$°(dn x R+). 

In the case of the wave equation the situation is more complicated. How- 
ever, it is known that given R w f\[o : T] f° r any / G C^°(dfl x [0,T]) we can 
construct R w f\-R + for all / G Co°(<9$7 x R + ) as soon as T > 2 max{d(x, dil) : 
x G $7} (see |KL2 |, [KKL|). We note that in this case it is also possible to 
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solve the inverse problem for the wave equation directly, without continu- 
ation of data. For this fact we refer to e.g. [BKa], | KKLfl , and, for a 
more general case of a wave equation which corresponds to a non-self adjoint 
symbol a(x,D), to [KL1|. 



5.2 Anisotropic boundary form and energy flux on a finite 
time interval. 

Let us consider the inverse problem for the wave equation with finite time 
observations. It is known that with appropriate boundary data measured on 
a time-interval [0, T] it is possible to reconstruct the operator in the domain 
{x : d(x,dQ) < T/2}. In the reconstruction procedure it is actually not 
necessary to know the Dirichlet-to- Neumann map R w f\\o l T] but only the 
antisymmetric form corresponding to it (see KKIfl ). This antisymmetric 
form, B T [f, h] is given by 

B T [f,h}= [ [ (f(x,t)R w h(x,t)-h(x,t)R w f(x,t)) dtdS x 
Jan Jo 

with real-valued f,h € C >°°(dn x R + ), where C°>°°(dn x R + ) consists of 
C^-functions with support in R + . As this antisymmetric form is related to 
the symplectic structure on L 2 (dil x [0, T]) it is often a more natural object 
in inverse problems than the Dirichlet-to-Neumann map. 

Next we consider the relations between B T and the finite time energy 
flux, 

U T [f,h] = J Re / / (R w f(x,t)d t h(x,t) + d t f(x,t)R w h(x,t)) dtdS x (42) 
* Jan Jo 

also considered on real- valued / and h. Hence, Tl T [f,f] = E w (ut \T) is 
the amount of energy that has passed into f2 through its boundary until 
t = T. In particular, we consider the question of finding the antisymmetric 
boundary form B T [f, h] from the energy fluxn T [/,/i] of form (EE) given for 
f,h& C '°°(dS} x R+). 

For f,h€ C°>°°{dn x R+), d t R w f = R w d t f so that 

B T [d t f, h]= I f T d t f(x, t)R w h(x, t) dtdS x + (43) 
Jan Jo 

+ I f d t h(x,t)R w f(x,t)dtdS x - [ h(x,T)R w f(x,T)dS x = 
Jan Jo Jan 

= 2U T [f, h]- f h(x, T)R w f(x, T) dS x . 
Jan 
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As any / G x R+) is of form / = d t f , f G C >°°(dtt x R + ), rep- 

resentation (|43|) makes it possible to evaluate B T [f, h] in terms of n T [/o, /i] 
as soon h(x, T) = 0. Similar considerations show that T1 T determines also 
B T [f, h] when f(x,T) = 0. We need also the following result proven in 
Appendix 2. 

Proposition 1 Let F,H£ C°°(dQ) be real-valued and a > 2. Then 

B T [t a F, t a H] = (U T [t a F, t a+l H] - U T [t a+l F, t a H}) . (44) 

Because any / £ C°'°°(dQ x R + ) can be represented as 
f(x, t) = f (x, t) + t a F(x), f (x, T) = 0, a > 2, 
and similar representation takes place for h, we see: 

Proposition 2 Knowing H T [f,h] for all f,h G C 0,oo (9f2 x R + ) we can 
determine B T [f, h] . 

5.3 Riemannian manifolds. 

We mention that all previous results can be directly generalized to man- 
ifolds. Indeed, let (M,g) be a Riemannian manifold with — A g being its 
Laplace-Beltrami operator. The corresponding Schrodinger operators on 
the manifold M are operators of the form A = — A g + q. Then the previ- 
ous arguments remain valid for Schrodinger operators on M. In particular, 
Theorem [l] claims that all data i.-vi. (and also vii. when zero is not an 
eigenvalue) are equivalent. 

5.4 Gauge transformations 

When we study equivalence of boundary data for general 2nd-order operators 
instead of the Schrodinger operators (Q)-(0), we have to take into account 
possible coordinate and gauge transformations. A general theory of gauge 
transformations in inverse problems for an arbitrary self-adjoint 2nd-order 
operator 

m m 

Au = ^2 a jk {x)djd k u + ^ b j (x)d j u + c(x)u, (45) 

j,k=l j=l 

v(A) = H 2 {n)nH^{n) 
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is developed in [KKL], see also [|K1[, [KK|. In this subsection we will con- 



sider a subclass of operators (|4j|) and show how the results obtained for 
Schrodinger operators can be generalized for them. Naturally, when study- 
ing a subclass of operators (fH|), we have to consider the subgroup of the 
complete group of the coordinate and gauge transformations which preserves 
this subclass. This subgroup is called the admissible group. 

As an example, we consider the class of the anisotropic conductivity 
operators, 

m 

a(x,D)v(x) = ~Y J d j (a jk (x)d k v(x)) + q(x)v(x) (46) 
j,k=l 

with the Dirichlet boundary condition (g), rather than a Schrodinger oper- 
ator (|l|). We note that all operators of form (^) that are self-adjoint with 
respect to the Lebesque measure on Q can be written in form (^). Then 

m 

{X h B(pt\ m : 1 = 1,2,...}, Btpi = a jk Ujd m , (47) 

j,k=i 

is the corresponding boundary spectral data. Similarly, we can define ellip- 
tic, hyperbolic, etc. boundary data of form ii.-vii. with the boundary oper- 
ator B of form (p7|). Although it seems that the question of the equivalence 
between data i.-vii. for operators fl46|) , ( f47| ) is completely analogous to that 
for Schrodinger operators, (|l]), (Q) there is a significant difference. Namely, 
it is, in principle, impossible to determine the corresponding Dirichlet-to- 
Neumann map, 

/ -> Bu[\q M 



with B of form (47) and also the corresponding non-stationary Dirichlet-to- 
Neumann maps R w , R h , and R s from the boundary spectral data ((!?]). As 
we will see below, the reason for this is that the admissible group of gauge 
and coordinate transformations which does not change the type of an oper- 
ator under consideration can change the boundary data. More precisely, let 
X : Vt — > X\qq = id\gQ be a diffeomorphism (coordinate transformation) 
in the domain S7. It gives rise to the transformation, Sx ■ L 2 (£l) — > L 2 (Q), 

Sxu(x) = u(X(x)). 

Similarly, any k G C°°(Q), k > in £1, gives rise to the gauge transformation, 
S K : L 2 (Q) L 2 (n), 

S K u(x) = k(x)u{x). 
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Then, S K o Sx '■ L 2 {0) — ► L 2 (f2) induces a transformation for operators 
which maps a(x, D) of form (|45| ) to the operator 

5(x, D)u = S K oS x a(x, D) (S^ o S^u) . (48) 

Next, assume that a(x,D) is of form (f5f). Taking 

k(X(x)) = |det(dX)| -1 (x), (49) 



we see that a(x,D) is of the same form ( [Iq ) as a(x,D) with, however, 
different coefficients 5^ 7^ ay, q ^ q. Furthermore, it can be shown that the 
only transformations of the form S K o Sx which preserve the form (46) of 



an operator in a given SI are given by ((4q)-(|49j), Moreover, if det(dX) = 1 
on dQ, the boundary spectral data of a(x, D) and a(x, D) coinside. On the 
other hand, in this case 

A z f = A z f + (BK)\ dn f, 

Analysing the proof of Lemma |l] in Appendix 1, we see that asymptotics 
( f2l|) remains valid for the operator ( |46| ) with, however, other factors a(x) 
and K(x) instead of p(x) and H{x). These factors depend only on aij{x). 
We can then follow the arguments of section 2, i. — * ii. to show that the 
boundary spectral data determine the Dirichlet-to-Robin map 

/ -> Bu f z \gn = Bu{ + au f z \ dn , z G C 

with an unknown function a £ C°°{dM). Motivated by this we say that 
operators A^ and A^ are gauge-equivalent if there is a function a such 
that f = A^/ + af. Similarly, two bilinear forms Ili 1 ^ and ni 2 '* are 
gauge-equivalent if there is a function a such that 

n«(/,/) = ni 2 ) (/,/)+ / a(x)\f(x)\ 2 ds x . 

Then, using similar considerations as in sections 2-4 we obtain: 
Proposition 3 Inverse problems i.-vi. for anisotropic conductivity opera- 



tors (46) are gauge- equivalent, in the sense that the data i. determines the 



gauge-equivalence class of any data i.-vi. and vise versa. 

We would note that this phenomenon does not occur in the Schrodinger 
case because its admmissible group of transformations consists of only changes 
of coordinates, that is, Sx- 

In connection to the inverse problem, we would mention that it is possible 
to find a(x, D) (upto a transformation from the group (|48|)) from either the 
boundary spectral data (f4~7|) or the Dirichlet-to-Robin map A z +a(x)I , z S C 
with unknown a (see e.g. fti), jKK|], flKKLf ). 
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6 Appendix 1 

In this appendix we will prove asymptotic expansion (pi]). To this end we 
introduce boundary normal coordinates (y,n), y = (y , • • • ,y m ~ 1 ), n > 
in a vicinity of d£l. Let x £ O. Denote by n = n(x) the distance (in 
the metric g) from x to the boundary dCl, n(x) = d g (x,dQ). When x 
is sufficiently close to <9f2 there is a unique point y = y(x) with n(x) = 
d g (x,y(x)). Introducing some (local) coordinates (y ,y m ~ l ) on 80, we 
use (y 1 (x), ■ • • , y m_1 (x), ra(x)) as (boundary normal) coordinates of x. The 
length element dl in these coordinates has a special form, 

dl 2 = dn 2 + g a p{y, n) dy a dy^ , a, /3 = 1, • • • , m — 1. 

Then 

n o) = 2 9 « ( l °sg(y, «)) ln=n ( 50 ) 

is the mean curvature of the surface dQ m = {x £ Q : n{x) = n }. In 
particular, 8Qq = d£l. 

Let a(x, D) be of the form 

a(x, D)v = — A g v + q(x)v. 

We denote by = — id y , D n = —id n , D s = —id s the derivatives in bound- 
ary normal coordinates. The crucial step in proving (f2l|) is the following 
lemma. 

Lemma 1 The Dirichlet-to- Neumann map A z of form (j^p, (Q) for 

a(x, D) has an asymptotic expansion 

oo 

Az ~ ^^Vk(y,D y )(-z)~( k ~ 1 ^ 2 , (51) 

k=0 

where Vk(y,D y ) are differential operators of order k. Expansion (51) means 
that 

i 

\\k - E^(y'^)(-^)" (fc " 1)/2 ii^ +1{ ^)^(^) ^ c ^ m 

k=0 

in the domain | arg(— z)\ < 5tt, < 5 < 1. In particular, 

H(y) P 
2 

1 

2 

7> 3 (z/, A/) = -\pd n g aP ( y ,Q)D a Dp + h{y,D y 



V (y, D y ) = -p, V x {y, D y ) = (52) 

V 2 (y,D y ) = -^pg^(y i 0)D a Df 3 +V 2 (y,D yJ , 
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where H(y) = H(y,0) and Vi, V3 are differential operators of order 1 with 
coefficients depending on the derivatives of the metric tensor and the poten- 
tial on dO,. 

Proof. By adding, if necessary, a positive constant to q{x) we can 
assume that o~(A) C {z : Re z > 0} where A the operator with symbol 
a(x,D) and the Dirichlet boundary condition. 

In [LU] it is shown that the Dirichlet-to-Neumann map for an elliptic 
differential operator of the second order is a pseudodifferential operator. 
Its symbol depends on the derivatives of the metric tensor at the bound- 
ary. To apply this approach, we introduce an auxiliary elliptic operator 
a(x, s, D X ,D S ) = a(x, D X ,D S ), 

a(x, D x , D s ) = a(x, D x ) — dg in (x, s) G f2 = $7 x R. 

Rewriting a(x, D x , D s ) in coordinates (y,s,n) we have 

a(x,D x ,D s ) = a(y,n,D y ,D n ,D s ) 

= -dl -d 2 s - g at3 {y,n)d a dp - H(y,n)d n + f a (y,n)d a + h(y,n). 



By [LU], a(x, s, D x , D s ) has a factorization 

a = (D n - iH(y, n) - iL(y, n, D y , D s )){D n + iL(y, n, D y , D s )), (53) 

The operator L(y,n, D y , D s ) is a pseudodifferential operator of the first order 
with respect to (y, s) which depends smoothly on the parameter n. The 
symbol of L(y, n, D y , D s ) has the form 

L(y,n,£,r) ~ ^ L m (y,n,£,r), y G R m ~\ n £ R+, 

where £ G R m_1 , r G R are dual to y and s, correspondingly. L m (y, s, £, r), 
m = 1, 0, — 1, . . . , are positive-homogeneous symbols of order m with respect 
to (£,t) and we take the principal symbol Li(y,n,^,r) to be negative. 

Using the calculus of pseudodifferential operators and comparing the 
terms of the same homogeneity with respect to (y, r) on the both sides of 
equality , we obtain 



^OjLjD^Lk + (d n + H)L m = a m , (54) 
j,k<i 7 ' 

j+k— \^\=m 
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where a m = a m (y, n, D y , D s ) is the term of homogeneity m with respect to 
(D y ,D s ) in a. In particular, a m = for m < 0. Using equation (|5^ ) we 
can find L m {y, n,£,r) in terms of coefficients of operator a. In particular, 
equation (^) for m = 2 implies that 



Ll(y,n,£,T) =-T} = -^1] 2 + T 2 , 1] = 

Further equations (|54|) give us recurrently that 



(55) 



L r 



1 

2^ 



rre+l<j,fc<l 
\j+fe-[7|=m+l 



(56) 



Formula fl56| ) implies that 



L m (y,n,£,T) = ^2 Pm-jiv^'Ov 3 > m<0, (57) 

jo(m)<j<m 

where Pj, (y, are polynomials of order fc with respect to £. Using the 
binomial Taylor expansion of 77 when r — > 00, 



m 

i=o 



i 



(58) 



where (?) = a(a — 1) . . . (a — fe + we see that 

L m (y,n,C,r) = ^p^ i (y,n,£)-r'. 



(59) 



Here pi (y, s, £) are again polynomials of order k with respect to £. These 
polynomials can be found by direct computations. In particular, from form 
( |55| ) for Li we see that 



(60) 



pfW,£) = -1, pL%,n,£) = ~^(y,n)Ufi- 
Furthermore, analysing formula (E(f) with m = we obtain that 



~(0), 

Po iv 



,n,0 = "^%^, P ( _%,n,0 = -7 dn^MUp- (61) 
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Since for L\ with the negative principal symbol, 

L(y,0,D y ,D s ) = Aq, 
where A z is the Dirichlet-to-Neumann map for a, i.e. 

A z / = 9n^|aQ, where a(x,D)vl = zv(, v(\oq = f, 



formulae (|57|), ( p9|) (with n = 0) determine the symbol of Ao- 

Now we are ready to derive the asymptotics of the Dirichlet-to-Neumann 
map A_ T 2 of the operator A when r -> oo, r G R. Consider the boundary 
value problem 



a(x,D x ,D s )v = 0, in £1, 



Its periodic solution v(x, s) = v(x)e lTS has the function v(x) being the solu- 
tion of the Dirichlet problem 

a(x, D)v + t 2 v = 0, in fi, 
v\dn = /• 

Since Bv = p §^| n=0 , we have 

A_ T .f = pe- lTS A (e lTS f). (62) 

Therefore, A_ T 2 is a pseudo differential operator with symbol 

A_ r2 (y,£) = p(y)L(y,0,£,T). (63) 

Summarizing the previous considerations, we obtain an asymptotic expan- 
sion for the operator A_ T 2, 

i 

A_ r2 = Y,n(y,D y ) T - k+1 + E l (y,D y ;T), (64) 

k=0 

where Vk{y,D y ) is a differential operator of order k and Ei(y, D y ;r) is a 
bounded operator from H l+1 (dQ) to L 2 (dQ) with the norm \ \Ei\\ < C\T~ l '. 
The exact form of the operators Vk(y,Dy) can be obtained by direct cal- 
culations using formulae (|55l)-([59|), (|62^), (p^). In particular, the first four 
polynomials have form (|52^). □ 

Lemma |l| yields immediately the following result. 
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Corollary 1 The derivative of Dirichlet-to- Neumann mapping, d z A z has 
the expansion 



1 , s_i 10 £kk—l 



d z A z ~ -pi-z)' 1 ' 2 + —V k {y,D y ){-z)-^ 2 , (65) 

k=2 

where Pk(y,D y ) are differential operators given in Lemma [J. In particular, 
the map d z A z determines the weight p and mean curvature H{y). 

Proof. From the principal symbol of d z A z we obtain p. It is then clear 
from formula ( p2|) for Vs(y,D) that the map d z A z determines the mean 
curvature H(y). □ 



7 Appendix 2 

In this appendix we give the proof of Proposition |l|. We can represent 
B T [t a F,t a H] in the form 

B T [t a F, t a H] = - [ T [ t a (8 u u taF H-F d u u taH ) dtdS x = 

— U f T [ (d u u ta+lp d t (t a H) + d t (t a+1 F)d u u taH ) dtdS x - 
a + 1 Jo J d n y J 

d u u ta+lF (x,T)H(x)dS x . (66) 

a + 1 Jan 

Similarly, 

B T [t a F, t a H] = — — [ T [ ( d v u ta+lH d t (t a F) + d t (t a+1 H)d u u taF ) dtdS x + 
a + 1 Jo Jan v ' 

+ ^7t/ d v u ta+lH (x,T)F(x)dS x . (67) 
a + 1 Jan 

Subtracting (|67|) from (|66|) and using definition (^g) of II T , we obtain that 

/ (d v u ta+lF {x, T) H(x) + d u u ta+lH {x, T) F{x)) dS x = 
Jan v ' 
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2T - a (n T [t a+1 if, t a H] + n T [t a F, t a+1 H]) . (68) 

On the other hand, denoting by r w {t, x,y), x,y G 50 the Schwartz kernel of 
the operator R w and using the fact that r w {-,x,y) = r w {-,y,x) we see that 

/ d u u ta+lF (x,T)H(x)dS x = [ d u u ta+lH (x,T)F(x)dS x . 
Jan J 80, 

Therefore, equation ( |68[ ) shows that 

I d„u ta+lF (x,T)H(x)dS x = 
Jan 

= T- a (n T [t a+1 F, t a H] + n T [t a F, t a+1 H}) . (69) 

Equation @ follows from (||) and (||). 
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